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EXECUTIVE SUMMARY 


Clusters of metal atoms are of Importance as models of substrates 
In studies of Chemisorption on metal surfaces and of catalysis of chemical 
reactions; they are also of Increasing Interest as novel molecular species 
In their own right that, surrounded by certain ligands, are possibly a whole 
new class of catalysts. Moreover, clusters of metal atoms are useful as 
prototype atomic-level host systems for the development and testing of models 
on Interatomic forces. These Interatomic forces are a basic Ingredient to 
the much more complex atomic modeling of the mechanical properties of metals 
that are used to address Important materials problems such as the effect of 
hydrogen on crack propagation. 

The approach we have taken to the study of metal atom clusters Is 
based on the application of ab Initio quantum mechanical approaches. This 
final report discusses our current research effort In the application of 
these methods. Because these large "molecular" systems pose special practical 
computational problems In the application of the quantum mechanical methods, 
there Is a special need to find simplifying techniques that do not compromise 
the reliability of the calculations. Our current research Is therefore 
directed towards various aspects of the Implementation of the Effective Core 
Potential (ECP) technique for the removal of the metal atom core electrons 
from the calculations. This final report discusses our recent progress In 
this 


area. 



QUANTUH CHEMICAL CALCULATION OF THE 
EQUILIBRIUM STRUCTURES OF SMALL 
METAL ATOM CLUSTERS 

by 

Luis R. Kahn 
INTRODUCTION 


The goal of this research Is to gain quantitative fundamental 
understanding about the properties of dusters of metal atoms. The In- 
creasing Interest In these novel "molecular" species stems from their use 
to model substrates In the study of chemisorption, and from the prospects 
that such studies will yield the long-sought understanding of the mecha- 
nisms of heterogeneous catalysis of chemical reactions. Moreover, these 
metal atom clusters surrounded by carbonyl groups, are Indeed being thought 
of as new types of catalysts In their own right. In every case, however, 
there Is still very little known yet about these Interesting new molecular 
species. Since the atomic level properties of these clusters are not yet 
readily available from experiment, the most reliable approach to their 
determination lies In the development and application of predictive ab Initio 
quantum mechanical approaches. 

This final report describes our current research effort designed 
to continue the development of a reliable, more accurate, and more efficient 
quantum mechanical method necessary for the calculation of the properties of 
clusters of metal atoms. The method deals with the Important simplification 
In the quantum mechanical calculation of metal atom clusters of the removal 
of the metal atom core electrons from the calculations. The method Is the 
Effective Core Potential (ECP) approach to simulate the effect of the core 
electrons on the valence electrons. 

SCIENTIFIC PROGRESS 


A New Theoretical Analysis of the ECP Method 

A long standing issue in the Approach of Effective Core Potentials 
(ECP's) has been the absence of a theoretical analysis on the basis of which 
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the properties of the ECP's could systematically be Improved* a theoretical 
analysis that would unify the evolving but seemingly disparate proposed 
improvements to the approach. We have succeeded in finding this needed unify 
ing theoretical analysis for the ECP approach, and a paper describing this 
work will be appearing in the International Journal of Quantum Chemistry. 

A preprint of this paper is found in Appendix A of this final report. 

We sunmarize here the background to this theoretical development. 

As is well known, the first ECP's were based on valence pseudo-orbitals con- 
structed from linear combinations of atomic orbitals of the same orbital 
angular momentum. While these first ECP's led to some excellent comparisons 
with molecular all electron calculations, they also led to some perplexing 
failures. An ad-hoc remedy was proposed by Hay et al,^^^ based on the dis- 
covery of the sensitivity of the molecular potential energy curves to the 
long-range properties of the ECP's. Whereas the ad-hoc remedy was found to 
lead to molecular potential energy curves that were in much better agreement 
with all electron results, a fundamental explanation was lacking. A funda- 
mental understanding of the problem was found in the recognition of the 
charge redistribution that was unwittingly being introduced in the valence 
pseudo-orbitals by the li near-combi nation-of -atomic-orbitals aspect of its 
definition. The improper long-range behavior of the ECP's was in retrospect 
only a symptom of this more fundamental deficiency in the valence pseudo- 
orbitals. The result has been the development of the "shape consistent" 
valence pseudo-orbitals. Whereas it appears that a number of workers dis- 
cussed this idea.^^"®^ Christiansen et al^®^ were^the first to demonstrate 
the important implications of "shape consistency" for the reliable calcula- 
tion of molecular potential energy curves. Moreover, the practical proce- 
dures of the Christiansen et al approach for the construction of "shape" 
consistent" valence pseudo-orbitals appear to be the most widely adopted 
presently. 

Notwithstanding the practical success of the "shape consistent" 
approach, a conmon concern has been not only that the core segment of the 
valence pseudo-orbital is arbitrary except for the single condition of 
normalization, but that no physical criteria had been established to decide 
even in principle just what additional conditions the core segment of the 
valence pseudo-orbital should satisfy. This concern has been heightened 
by the observation in practice of unexpectedly sharp spatial characteristics 
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In the "shape consistent" ECP's. In the absence of a comprehensive theoret- 
ical analysis It has been difficult to decide whether these properties of 
the ECP are physically Important properties or merely artifacts of the In- 
sufficient physical definition of the core segment of the valence pseudo- 
robltal. 

Through the valence pseudo-orbitals and the ECP one Is, of course, 
attempting to simulate as closely as possible the parent atomic all-electron 
Hartree-Fock solution. Indeed, Rappe et al^^^ have proposed that the core 
segment of the valence pseudo-orbital also be required to minimize the error 
In the valence-valence Interactions. Their approach, named the "hamlltonlan 
and shape consi-tent" approach, points to an Important additional physical 
criterion on the basis of which to remove some of the arbitrariness left In 
the definition of the core segment of the valence pseudo-orbitals. The 
minimization of the valence-valence Interactions Is a global property of the 
orbitals and, notwithstanding its importance. Is a quite disparate property 
of the parent all electron Hartree-Fock orbitals than the normalization condi 
tion. Moreover, the details of the valence-valence interactions are closely 
dependent on the particular electronic state, and change from one atom to the 
next In a complicated manner. 

Our solution to these Important Issues In the ECP approach is de- 
scribed in detail In Appendix A. We point out here the main properties and 

f21 

conclusions of this novel development.' ' 

(1) A new set of attributes of the atomic orbitals Is 
found that alone determine the valence-valence Inter- 
actions fully. These new physical measures of orbital 
characteristics are found In the "moment accumulation" 
functions defined as 

O'” 

where ”nt and are atomic radial functions. 

(2) Since these "moment accumulation" functions determine 
the valence-valence Interactions, the minimization of 
the error In these orbital properties In turn mini- 
mizes the error in the valence-valence interactions. 
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The latter is shown to be equivalent to the require- 
ment that the core segment of the the valence pseudo 
orbital satisfy the conditions 



The point at which the core and valence segments of 
the valence pseudo-orbital are joined is denoted as rj^ 

In the above, the superscript tilde indicates the 
terms constructed from valence pseudo-orbitals 
rather than the parent valence Hartree-Fock orbitals. 

The set of valence orbitals is denoted as V. 

(3) The previous sole requirement of normalization is found 

now to be only one of the new conditions. The new requirements, 
however, are also shown to be discrete individual condition* on 
the orbitals of the same succint nature as the 
normalization condition. Since the satisfaction 
of these conditions minimizes the error in the 
valence-valence interactions, the "Hamiltonian 
and shape consistent" condition is, through these 
conditions, shown in fact not to be a disparate 
set of conditions from the normalization condition. 

Moreover, the satisfaction of the "Hamiltonian and 
shape consistent" condition is hereby extricated 
from the complexities of the state dependent atomic 
energy expressions. 

(4) Indeed, it is possible to extend the required set 
of conditions systematically to attain an ever in- 
creasing simulation of the orbital characteristics 
(that, in turn, determine the interactions among 
electrons) without being limited by the difficulties 
of finding some pertinent model energy expression. 

There is, of course, a limit to such a generaliza- 
tion imposed by our overriding goal to project out 
core states through the requirement of nodelessness 
in the valence pseudo-orbital, as well as our prac- 
tical requirement of obtaining valence pseudo-orbitals 
that are "smooth" over the core region. 
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(5) A novel expression for the ECP Is obtained In terms 
of the "moment accumulation" functions that allows 
us to prove that the residual In the valence-valence 
Interactions strictly vanishes for r ^ rj||| If the 
valence pseudo-orbital satisfies the "hamlltonlan 
and shape consistent" conditions. This explains In 
large measure the extremely sharp behavior already 
observed in "shape consistent" ECP's around 
Whereas the residual In the valence-valence Inter- 
actions Is present for < »*|nt the residual Is forced 
to vanish for r ^ rjj^ by virtue of the conditions 
built Into the core-segment of the valence pseudo- 
orbital. This effect takes place at whatever value 
rj^ is given, and, therefore, explains the observed 
change In the shape of the ECP's with changes of r^^. 

In summary, this new theoretical analysis yields understanding of 
previously observed yet unexplained properties of the "shape consistent" 
ECP's, it unifies the previous "shape consistent" approach with the more 
recent "shape and hamiltonian consistent" approach, it shows that the orbital 
attributes that control the valence energy are the "moment accumulation" 
functions, it affords an analysis of the long-range behavior of the ECP in 
terms of the "moment accumulation" function, and it provides a foundation 
for the future construction of systematically improved ECP's. 

Evaluation of Matrix Elements of the ECP 


In the semiannual progress report of October 23, 1981, we presented 
the reformulation of the integrals over cartesian gaussian basis functions 
(GTO) of the molecular ECP so that the shell structure of the GTO's could be 
exploited. A computer program that yields all the terms In these formulae 
shell by shell has been prepared and a copy is given in Appendix B. This pro- 
gram yields all the unique non-zero terms, and Is thus the replacement for the 
ENCODER routines of the current SPDF ECP integral program. The other major in- 
put to these formulae are the radial integrals, I^ ^ and The routines 

that calculate these latter integrals are the same ones that are already 
present in the current SPDF ECP Integral program. 
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Whereas the calculation of the ECP integrals over shells Is likely 
to yield a large computational savings, there are other efforts under way by 
workers at the National Bureau of Standards that have promise to Improve 
the efficiency of the evaluation of the radial ^ and Integrals. 

This work alms at the direct numerical calculation of the Integrals over 
the ECP on a grid of points. This approach may find advantages In the 1 'iter- 
change of the sum over terms with the sum over quadrature points. Another 
goal of the latter work Is to bypass the step of fitting the numerical ECP 
with GTO's. This step has become one of Increased difficulty as a result 
of the very sharp approach to zero of the "shape consistent" ECP's as r 
approaches r^^; this property of the ECP can be expected to remain In the 
"hamiltonian and shape consistent" approach on the basis of the analysis 
presented above In terms of the "moment accumulation" functions. 

There remain, thus, well defined technical problems In the Imple- 
mentation 0 ^ various Ideas for the Improved efficiency of the ECP Integral 
calculation. The prospects are gooo, however, that these will be solved, 
and that thereby the use of the ECP approach will be extended from a research 
tool to a bonafide "engineering" tool for chemical research. 
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ABSTRAQ 

A one-component approach to molecular electronic structure Is 
discussed that Includes the dominant relativistic effects on valence 
electrons and yet allows the use of the traditional quantum chemistry tech- 
niques. The approach starts with one-component Cowan-Orlffin relativistic 
orbitals that successfully Incorporate the effects of the mass-velocity and 
Darwin terms present In more complicated wavefunctlons such as the Dlrac- 
Hartree-Fock. The approach then constructs "relativistic" effective core 
potentials (RECP's) from these orbitals, and uses these to bring the 
relativistic effects Into the molecular electronic calculations. The use 
of effective 1-eIectron spin-orbit operators In conjunction with these one- 
component wavefunctlons to Include the effects of spin-orbit coupling Is 
discussed. Applications to molecular systems involving heavy atoms and 
comparisons with available spectroscopic data on molecular geometries and 
excitation energies Is presented. Finally, a new approach to ’'.he construc- 
tion of RECP's encompassing the Hamiltonian and shape consistent approach 
Is presented together with a novel analysis of the long-range behavior of 
the RECP's. 


1. INTRODUCTION 


This paper reviews various aspects of a practical but nontraditlonal 
approach to molecular relativistic electronic structure calculations. Mhere- 
as It Is traditional to approach relativistic electronic structure In terms 
of four-component^^'^^, or even two-component^^*^^* wavefunctlons, the 
present approach Is based on one-component wavefunctlons familiar from the 
nonrelativistic theory. It appears* moreover, that the relativistic contri- 
butions to molecular properties of chemical Interest are being reliably calcu- 
lated In all cases examined by this approach despite the approximations from 
which the simplicity of the approach derives. 

2. APPROACH 

This one-component approach to molecular relativistic electronic 
structure calculations has two basic Ingredients. These are: 

(1) One-component Cowan-Griffin relativistic 
atomic orbltals^®^ 

(2) "Relativistic" effective core potentials 

We review the essential properties of these two Ingredients to the present 
approach next. 


2.1 Relativistic Atomic Orbitals 

The one-component relativistic atomic orbitals that are the basis 
of this approach are obtained from the Cowan-Griffin equations. The 
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latter may be derived qualitatively from the local potential approximation 
to the D1rac*Hartree-Fock equations 


+1|8 -VW+±1Q^ 0) 

+ -J-a* (2) 


where P,^ and Q,^ are the "large" and "small" components, respectively, 
where a » 1/137.036 Is the fine-structure constant, and where 


H-t If i**"t 

K«-t-l If + t 


( 3 ) 

(«) 


Solving Equation (1) for Q,^ and Qj^, inserting these In 
Equation (2), and averaging the explicitly j-dependent K term. 



K 


ar - I 


( 5 ) 


ORIGINAL A * ' 

^ OF POOR QUALU Y 

one obtains an equation that In Its form Is the Cowan-6r1ff1n equation. 
The actual CoMan-6r1ff1n, or relativistic Hartree-Fock (RHF), equation Is 



( 6 ) 


where contains the non-re1at1v1st1c terms found In the Hartree-Fock 
equation. 



' L 

l it* 2r» 



( 7 ) 


whereas and Hp are the relativistic terms known as the mass-velocity 
and Darwin terms, respectively. 



The potential In Equation (7) contains the usual nonlocal Hartree-Fock 
potential. However, the analogous potential In Equations (8) and (9) Is 
replaced by a "local exchange" potential, 
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The atomic orbitals obtained from the RHF equation show the 
typical qualitative contraction of Iom i quantum nund>er orbitals and the 
expansion of the high i quantum nim^er orbitals when compared to the non- 
relativistic orbitals. This Is Illustrated In Figures 1 and 2 where the 
comparisons are for the 7s and 5f orbitals of Uranium atoms, respectively. 
One also finds good quantitative agreement between the orbital properties 
obtained from the average of Dirac-Hartree-Fock (DHF) and RHF orbitals. 

This Is Illustrated In Table 1 for both orbital radial characteristics and 
orbital energies In the Uranium atom. Moreover, virtually Identical 
excitation energies are shown In Table 2 for the DHF and RHF calculations 
on the Au atom, whereas the nonrelativistic results predict even the 
Incorrect ordering of the states for the 5d -► 6s and 6s -*• 6p transitions. 

2.2. "Relativistic** Effective Core Potentials 

The second Ingredient to the approach Is the use of Effective 
Core Potentials (ECP*s) to bring the relatl^^lstlc effects on the valence 
electrons from the atomic calculations Into the molecular calculations. The 
ECP approach was first developed to serve as a device to bring Into a 
molecular calculation the effects onto the valence electrons of the 
chemically Inert core electrons. Since the sources of the relativistic 
direct and Indirect effects on the valence electrons also are localized 
about the atomic nucleus, the role of the ECP*s was expanded to bring the 
relativistic effects from the atomic calculations Into the molecular calcu- 
lations as well. In this expanded role the ECP became an "relativistic" 

ECP or RECP.^^^ 
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FIGURE 2. COMPARISON OF THE NON-RELATIVISTIC (NR) AND RELATIVISTIC (R) 
HARTREE-FOCK 5f ORBITALS IN URANIUM ATOM 
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Table 1. Coaparlaon of Orbital Sixat and Cnergiaa for tba Uraniua 
Atoa. (Raf . 6 ) 



HF* 

rhf'* 

DHF® 

imp 

(average) 



<T*> (bohr*) 



5f 

5f; 

1.94 

2.57 

2.53 

2.67 

2.61 

6d4. 

9.63 

12.4 

11.1 

13.7 

12.7 

7a 

28.8 

21.5 

21.8 

21.8 

®P- 

bp+ 

4.11 

3.73 

3.13 

4.08 

3.76 



Orbital Energy (a.u.) 



5f 

5f; 

-0.634 

-0.331 

-0.352 

-0.297 

-0.320 

6d 

6d; 

-0.267 

-0.188 

-0.208 

-0.172 

-0.186 

7a 

-0.167 

-0.201 

-0.199 

-0.199 


-1.04 

-1.086 

-1.363 

-0.959 

-1.094 


^onrelativlatlc Hartree-Fock. 

Relativistic Hartree-Fock (Cowan-Griff In). 
®Dlrac-Hartree-Fock. 
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Table 2. Excitation Energies for the Au Atom. (Ref. 9) 


Excltetloo Energy (cV) 

W MF 




ORIG’ 
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We review the basic properties of the RECP's next, and for this 
purpose confine the Illustrative example to a nonrelativistic case. Con- 
sider for the sake of Illustration an atom with a single valence electron 
such as the LI atom. The atomic equation for the valence orbital Is written 
schematically as 


I • w*') 

I tr* 



U ) P - 6 P 


(10) 


where G physically represents all the Interactions with all other electrons 
which. In this Illustration, are the core electrons. This equation Is 
assumed to have been solved. We define a local potential to simulate the 
effect of the nonlocal U by^^®^ 


..locm 

U ir) 


p 


(11) 


with the property that It ensures by construction that, given D and P, 
Equation (10) Is simulated by 


I 1 JL 4. iM.l 

ri ar^ ^ 



( 12 ) 


We show In Figure 3 a nodeless valence orbital (P/r) corresponding 
to the ?s orbital In LI atom.^^®^ This figure also shows another choice 


AMPLITUDE (A.U.) 
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FIGURE 3. A MODELESS 2s ORBITAL IN LITHIUM ATOM. AND THE CORRE- 
SPONDING NOOEFUL 2s HARTREE-FOCK ORBITAL^s 
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for the 2s orbital which contains a node. We want to gloss over the 
Important con^llcatlons of the various other choices of P here and focus 
on the nodeless one In order to maintain the s1n^)11c1tv of the discussion. 
Using the nodeless 2s orbital one obtains the potential shown labeled as 
In Figure Indeed a different potential Is obtained depending on 

whether the valence orbital Is the 2p or 3d orbital, etc., and this yields 
the other potentials In Figure 4 labeled as Up and It Is found 

that each such potential Is the appropriate one for the spectrum of states 
of Its own angular momentum as Is schematically Illustrated In Figure 5. 
This points out the Important angular momentum dependence of the local 
potentials, and It explains the construction of the net ECP In terms of 
these local potentials confined with angular momentum projection operators 
« tn(”> 


grP • IPC|L 

u -II IU>U cr)<t^| 

In tM< ' 


(13) 


The effects onto the valence electrons embodied In the ECP's are 
brought Into the molecular calculations via the effective valence electron 
molecular hamlltonlan. 





(lA) 


Potenfiol (Horfree) 
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FIGURE 4. THE s. p, AND d LOCAL EFFECTIVE CORE POTENTIALS FOR THE 
LITHIUM ATOM 



FIGURE 5. SCHEMATIC REPRESENTATION OF THE ONE-ELECTRON ORBITAL SPECTRUM 
ASSOCIATED WITH THE s. p, and d LOCAL EFFECTIVE CORE POTENTIALS 
FOR THE LITHIUM ATOM, RESPECTIVELY 
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Apart from the ECP, all other operators In Equation (14) are the usual non- 
relativistic operators, and they range only over the valence electrons. 

Thus, the only new element In carrying out a molecular calculation with 
ECP's Is the evaluation of Integrals of the ECP In a molecular basis. Com- 
puter programs to calculate these new Integrals routinely are found Inte- 
grated In various quantum chemistry packages. Whereas In special cases 
(e.g., diatomic molecules) It Is efficient to evaluate the Integrals using 
the ECP In Its numerical form^^^^, for polyatomic calculations In a gausslan 
basis this direct approach has as yet been found prohibitively Inefficient. 
It Is found, however, that If the ECP Is first fit with analytical gausslan 
forms such as 


Z «*t 1 -tj 1 , 05) 

(121 

then the multicenter Integrals can be efficiently calculated.' ' Therefore, 
polyatomic calculations with ECP's presently use the above noted expansions 
of the ECP's. 

A direct comparison of the valence molecular orbitals resulting 
from all-electron (AE) and valence-electron (VE) calculations for the 
molecules L1H and Ll 2 ^ was possible since the AE and VE orbitals 
could be calculated from directly comparable generalized Valence Bond wave- 
functions^ This con^arlson Is shown In Figures 6 and 7 for the LiH and 
Ll 2 ^ molecules, respectively. The comparison shows the AE and VE molecular 
orbitals to be virtually Identical thereby establishing confidence In the 
validity of the basic ECP approach. ^ 


Amplitude (A.U.) 




Distonce (Bohr) 


FIGURE 6. COMPARISON OF A VALENCE MOLECULAR ORBITAL FOR L1H X'z OBTAINED 
FROM THE VALENCE ELECTRON (VE) ECP CALCUALTION AND A COMPARABLE 
ALL-ELECTRON (AE) CALCULATION 
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FIGURE 7. COMPARISON OF THE VALENCE MOLECULAR ORBITAL FOR LI? xh 

OBTAINED FROM THE VALENCE ELECTRON (VE) ECP CALCULATION AND 
A COMPARABLE ALL-ELECTRON (AE) CALCULATION 



OKIC.. - i 
OF POOR C. 

17 

The relativistic ECP's, or RECP's, are based on the one-component 
Cowan-Griffin (RHF) equation. The RHF equation for a valence orbital P^^^ 
Is assumed to have been solved, and to have the form 



CORF RFL 

where u consists of the Interactions with the core electrons, U con- 

VAL 

sists of the relativistic mass-velocity and Darwin terms, and U consists 
of the Interactions with the other valence electrons. To construct the RECP 
from this starting point, one seeks to obtain a new valence orbital (pseudo- 
orbital), and a local potential, Uj^(r), that combine to satisfy the 
equation 


l.ii 

Xr*' 




( 1 ?) 


The local potential is defined as^^’*^^ 


U/rl . t - 
t 



(18) 


The U term differs from the term of Equation (16) in that, where 
appropriate, the valence interactions are constructed from the valence pseudo- 
orbitals, rather than the original valence orbitals, P^^j. 
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The detailed properties of the valence pseudo-orbitals, In 
Equation (18) are of critical Importance In determining the reliability of 
the RECP's. The RECP Itself may be thought of as a device for the embedding 
of the valence orbital properties. He shall discuss some of these orbital 
properties In greater detail In a subsequent section. Here we only point 
out three broad desired characteristics In constructing the valence pseudo- 
orbltals^^^’^^^: 

(1) Nodelessness 

(2) Maximum similarity with the parent valence orbital, 

(3) Minimal special characteristics In the core region. 

A schematic comparison of the parent valence orbital, P^^, with two possible 

-v. 

choices for the valence pseudo-orbitals, P^^^, Is shown In Figure 8. 

3. SPIN-ORBIT COUPLING 

The one-component approach to molecular relativistic electronic 
structure calculations reviewed here does not Include the effects of the 
spin-orbit coupling operator. The Important effects of spin-orbit coupling 
have to be obtained In a subsequent calculation In which the energies Ej 
and states Vj obtained from the RECP calculations are used to set up a net 
hamlltonlan matrix 



* 1,1 


<^li IV 


(19) 
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The hamlltonlan matrix Is then diagonalized to obtain the molecular energy 

A 

In the presence of the spin-orbit operator, V^q. 

Three levels of decreasing sophistication but also of 
decreasing computational difficulty have emerged In recent years for the 
calculation of the spin-orbit matrix ®l««ents,^'rj| Vjq|v^ These are; 

(1) Use of the rigorous spin-orbit Interaction of the 
Brelt-Paull hamlltonlan requiring the computation 
of one-electron and tMO-electron multi-center 
Integrals 

(2) Use of an effective one-electron operator of the 
form 



where is sum Is over the nuclear centers, A . 

(3) Use of the "atoms-ln-molecules" method. 

3.1. Atoms- In-Molecules 

The computationally simplest and quickest approach to obtain some 
semi quantitative measure of the effect of spin-orbit coupling on the 
molecular states calculated using the RECP's Is the atoms-ln-molecules 
approach. The approximation Is valid when the molecular states retain 
the same spin-orbit coupling as the parent atomic states. The spin-orbit 
matrix element^tj hsoh are assumed to depend on the spin-orbit splittings 
of the constituent fragments, and remain constant as the atoms are brought 
together to form the molecule. 
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To Illustrate this approach we consider the XeF molecule. 

The potential energy curves calculated with an RECP but prior to Including 
sp1n*orb1t coupling are shown In Figure 9. The atomic Ionic limit splits 
upon Introducing spin-orbit coupling 



♦ F'CS,) 

( 21 ) 

Xt (\1 + fi\) 


Let X be half of the spin-orbit splitting between the Xe^C^Py^) 
states. The atoms- In-molecules approximation yields the net hamlltonlan 
matrix, 



The molecular states obtained by diagonalizing this matrix are 
shown In Figure 10. In turn, the emission wavelengths to the ground states 
calculated from these curves are given In Table This table shows 

that the correspondence between the calculated and the experimentally ob- 
served emission wavelengths Is unambiguous. The ready assignment of the 
emission wavelengths that follows In this example shows the usefulness of 


I 2 3 4 5 6 

r{&) 

FIGURE 9. COMPARISON OF THE NONRELATIVISTIC (NR) AND RELATIVISTIC (R) 
VALENCE ELECTRON (VE) POTENTIAL ENERGY CURVES FOR THE LOW- 
LYING ELECTRONIC STATES OF THE XeF MOLECULE PRIOR TO THE 
INCLUSION OF THE SPIN-ORBIT EFFECTS 





FIGURE 10. THE RELATIVISTIC (R) VALENCE ELECTRON (VE) POTENTIAL ENERGY 

CURVES FOR THE LOU-LYING ELECTRONIC STATES OF THE XeF NOLECULE 
INCLUDING SPIN-ORBIT EFFECTS (USING THE ATOMS-IN-MOLECULES 
APPROXIMATION) 
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Table 3. Emission Wavelengths (nanometer) In XeF Molecule. ‘ 


TRANSITION 

, CALCULATED 

EXPERIMENT 


390 

352 


299 

269 


922 

950 


307 


zn - 1 

368 


«W.R. WADT. PJ. HAY. 

L.R. KAHN. J. 

CHER. PNYS. 6L 


1752 (1978) 



this simple approximation In giving a quick broad guide to the effects of 
spin-orbit coupling. Note that the emission wavelengths prior to Including 
the atoms-ln-molecules spin-orbit effects* also given In Table 3, provide no 
correlation with the experimental emission wavelengths. 

3.2. Effective One-Electron Spin-Orbit Operators 

The most rigorous approach to the Inclusion of spin-orbit effects 
Is the use of the full molecular spin-orbit term In the Brelt-Paull hamll- 



Thls Is also the computationally most difficult approach since It requires 

the computation of multicenter Integrals of both one-electron and two-electron 

M 9-221 

Integrals. Although such molecular computations are now feasible' . 
they remain the exception rather than the norm because of the associated 
computational difficulty. There Is therefore a need to find reliable 
approximations to Equation (23) that go beyond the atoms-ln-molecules approach 
In dealing with the changes that occur as the atoms are brought together* 
and that yet remain computationally practical. 

The Increasingly large numbers of chemically Inert core electrons 
as one moves down the periodic table are an even larger computational obstacle 
to the use of the rigorous molecular spin-orbit Interaction* Equation (23)* 
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than thay are for the calculation of the nonrelativistic molecular energy. 

The problem lies In the presence of the two-electron spin-orbit operators. 

There 1s» however » good reason to think that the dominant contribution of 
these two-electron spin-orbit operators Is contained In the Interaction 
between core and valence electrons mediated by this operator. Indeed, 
examination of values of all Integrals needed for the evaluation of the 
matrix element C^g| ^so| ^ig^for the okygen molecule, for example, shows 
the core-valence two-electron spin-orbit Integrals to be one order o f 

(21 1 

magnitude larger than the valence-valence two-electron spin-orbit Integrals.' 

It Is possible to Isolate the specific form of the two-electron 

core-valence spin-orbit Interaction. We take as example the manner In which 

(23) 

the analogous problem has been solved In the nonrelativistic case.' ' It 
is well known that the two-electron core-valence nonrelativistic Interactions 
can be expressed as a net potential due to the core electrons. 

hh-m ^ 

e 


acting on the valence electrons (the sum In Equation (24)) Is over all core 
orbitals). This potential Is combined with remaining valence operators to 
yield the nonrelativistic effective valence-electron hamlltonlan. 



+ 


i t 


I 
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where the operator Indexes range only over the valence electrons 

Applying the same approach to the rigorous spin-orbit Interaction. 
Equation (23). we start by partitioning the mahy-electron wavefunctlon. We 
divide the all -electron one-component wavefunctlon Into a Cowan-Griffin set 
of core orbitals. and an unspecified valence wavefunctlon orthogonal 

. . a«» 

to these core-orbitals. 9 t 




(26) 


where ^ stands for the antlsyirmetrlzer. Integrating the rigorous spin- 

orbit Interaction. C^q. over the core-electrons and combining the resulting 

core-valence terms with the valence-valence terms, we find the following 

(241 

effective valence- electron spin-orbit Interaction:' ' 



where the operator Indexes range only over the valence electrons. The 
remaining core-core terms are all zero. 
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The core electrons are showm by Equation (27) to formally 
contribute to the effective valence-electron spin-orbit Interaction via 
the gradient of the potential 

(28) 

Since this potential Is of opposite sign to the nuclear potential, the 
core electrons can be said to "shield" the nucleus In the spin-orbit Inter- 
action analogously as In the nonrelativistic energy. On the basis of the 
localized atomic nature of the core orbitals, the potential In Equation (28) 
may be expressed as a superposition of atomic potentials each centered on 
one of the nuclei. Combining these atomic potentials with the nucKar 
potentials to obtain a net potential on each nucleus, v^^^, one may reexpress 
Equation (27) as 



where again the operator Indexes range only over the valence electrons. 

The bare nuclei one-electron spin-orbit terms are found In practice 
to be two orders of magnitude larger than the valence- valence two-electron 
spin-orbit Interactions. The effective one-electron spin-orbit terms In 
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V$^* ^Quatlon (29) in turn contain the “shielded” bare nuclei terms. These 
op*'e1ectron terms are expected to be thus at least an order of magnitude 
larger than the two-electron spin-orbit Interaction remaln'ing In Equation (29). 
On this basis* It Is reasonable to expect approximate one-electron spin-orbit 
operators to account already for most of the chemically Important molecular 
properties. 

Whereas one-electron operators of some sophistication have been 
described recently In the literature, we wish to Illustrate here how 
well even simpler one-electron approximation appear to work.^^^^ This 
simplest of approximations Is based on letting 



(30) 


to obtain the following approximation to Equation (29), 





u 




(31) 


eff 

Z Is a parameter adjusted to match either theoretical or experimental 
atomic spin-orbit parameters. Moreover, In molecular calculations all but 
the one-center Integrals of Equation (31) are neglected. The latter con- 
ditlon can probably be removed by allowing the parameter Z to have a 
radial dependence reflecting the Increased shielding at larger distances 
from the nucleus. Table 4 shows spin-orbit parameters calculated for states 


f 
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Table 4. Comparison of Silicon atom spin-orbit constants (In cm*^) 
calculated by Stevens and Krauss (SK) (ref. 26) and by the 
effective one-electron one-center spin-orbit operator, 
(ref. 27) 


State 

VE 

(5s7p) 

SK 

SK(3+)* 

Expt 

SI [3s^3p^] 

148.9 

157.8 

152.7 

148.9 

Sl"’[3s^3p’] 

179.8 

201.4 

196.3 

191.3 

Sl*[3s^4p'] ^P 

32.7 

36,9 

36.6 

40.0 

Sl^^[3p^) ^P 

250.2 

— 

310.1 

306.9 


*Same as SK except that the spin-orbit operator Is derived from 
Si(3+) with the (ls*2^2pr) core orbitals frozen as In SI (^P). 
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of silicon atom Mherein the valence 3p orbitals undergoes large changes. 

A reasonable comparison Is still shown In Table 4 with the more sophisti- 
cated one-electron approximation as well as with experiment for all states 
(except SI ), notwithstanding the extreme simplicity of this effective 
one-electron approximation. Surprisingly good agreement with spin-orbit 
coupling matrix elements using the full Brelt-Paull hamlltonlan between 
electronic states of the rare gas oxides ArO, KrO. and XeO as a function of 
internuclear distance Is obtained by Langhoff^^^^ and shown In Figure 11. 

4. REPRESENTATIVE APPLICATIONS 


We briefly review next some representative molecular applications 
of the one-component approach to relativistic structure calculations dis- 
cussed here. 

To lllustrete the relativistic effect on chemical bonds, we 

compare relativistic (R) and nonrelativistic (NR) calculations of the 

4 (9) 

potential energy curves of the AuH and AuCT In their X'z ground states.^ ' 

The potential energy curves are compared In Figure 12 and Tables 5 and 6. 

It Is observed that the equilibrium bond length Is predicted to be smaller 

by the relativistic calculations In both cases. In AuH, for example, a bond 

contraction of 0.3 A Is predicted that yields a bond length differing by only 

0.01 A from the experimental value of 1.52 A. 

Since It has been traditional to correlate binding characteristics 

(gl 

directly with orbital properties, we have Interpreted' ' this bond contrac- 
tion as originating In the relativistic contraction of the Au 6s bonding 



FIGURE 11. SPIN-ORBIT COUPLING MATRIX ELEMENTS BETWEEN ELECTRONIC STATES 


or THE RARE GAS OXIDES COMPUTED BY LANGHOFF (REFERENCE 22) 
USING THE FULL BREIT-PAULI HAMILTONIAN (SOLID LINES) AND THE 
1 -ELECTRON 1 -CENTER EFFECTIVE SPIN-ORBIT OPERATOR (CIRCLES 
AND CROSSES) 
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R(A) 


FIGURE 12. THE POTENTIAL ENERGY CURVES FOR THE GROUND STATES OF AuH AND 
AuCI OBTAINED USING A NONRELATIVISTIC (NR) AND RELATIVISTIC 
(R) ECP AND A GENERALIZED VALENCE BOND WAVEFUNCTION 
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Table 5. Spectroscopic Properties of AuH molecule* 


Ad ix'l*) 

»,(A) 



Hoard. ECP 

Hr 

1.763 

0.99 

1387 

CVB-l 

1.820 

1.52 

1203 

roL>ci 

1.807 

1.57 

12.17 

lol BCP 

HP 

1.508 

1.55 

2014 

CVl'l 

1.514 

2.14 

1891 

P01/-C1 

1.522 

2.23 

1871 

Expel 

1.5237 

3.37 

2305 

1-CeotcT BlroC'Fock 

Noortl. 

1.745 

— 

2296 

Rcl. 

1.659 


2178 

«P. J. Hay. W. R. Wadt, L. 

R. Kahn, 

F. W. Bobrowicz , J. 

Che«. Phys., 


64, 984 (1978). 
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Table 6. Spectroscopic Properties of AuCl Molecule* 


AuCl 


D^(eV) 

«.(»-*) 

Nonrcl. ECP 
CVB-1 

2.447 

2.58 

277 

lei. ECP 

6VB-1 

2.283 

1.96 

298 

POL-a 

2.291 

2.39 

306 

Exptl 

— 

3.5 t 0.1 

382 

*P. J. Hey, H. R. 

Uadtp L. Re Kahn, 

F. H. Bobrowicz , J. 

Cheu. 


Phy*., M, 984 (1978) 
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orbital (<^>6s Is reduced from 3.7 Bohr to 3.0 Bohr). Ziegler, SnIJders 
and Baerends^^^^, hoMever, have established that the bond contraction first 
has Its origin In the relativistic mass velocity and Darwin terms prior to 
their effect on the Au 6d orbital. It Is Indeed reasonable. In retrospect, 
that the same relativistic terms that would further attract the Au 6s orbital 
to the Au nucleus, and cause Its relativistic contraction, should also 
attract the approaching hydrogen Is orbital, the other bonding orbital, and 
thereby lead to the bond contraction. 

The decrease In the binding energy In AuCI, an Ionic molecule. Is 
Interpreted to be a result of the deferral to smaller Internuclear distances 

+ - ( 9 ) 

of the onset of the admixture of the Ionic configuration Au Cl . The 
delay of this onset Is the result of the Increased energy of the Au^ + Cl’ 
asymptote relative to the covalent neutral atoms asymptote; It Is caused 
by the relativistic energy stabilization of the Au 6s orbital. 

Some of the Important and exciting research problems for the 
appl icablon of relativistic molecular methods lie In the chemistry of 
transition-metal complexes. Ue Illustrate the applicability of the 
present one-component approach to these types of molecules with some results 
by Hay^^®^ on the PtCl 3 (C 2 H^)’ complex (Zelse's salt). Figure 13 shows good 
agreement In a comparison of the calculated geometric parameters of Zelse's 
salt computed using a Pt RECP, and the geometric parameters obtained from 
neutron diffraction experiments. Moreover, the perpendicular orlentaton of 
the ethylene ligand relative to the PtClj’ plane Is correctly predicted to 
be the stable form with the coplanar configuration lying 15 kcal/mol higher 
In energy. Finally, the bending of the CH 2 groups away from the Pt atom Is 
also predicted by these calculations. 



37 


Cl 

/ 


Cl 


m 

Pi 


Cl 


/ 


H 

iee.4* 

y (ies.7») 


OF POOH QuH!..> i ^ 


Ml 

( 2 . 02 ) 


( 1.274) 


K 

A A • • 


0.0 kcQl 


Cl 


Cl 


Pf 


6 \ 


Ptcii 


V I76.9* 


2.42 




i' 


6 


ISO* 


L330 

(1.327) 


t 


4 14.9 kcol 


4 28.5 kcol 


FIGURE 13. GEOMETRICAL AND ENERGETIC PARAMETERS CALCULATED BY HAY (REFERENCE 
UfaUKt wu COMPARISON WITH EXPERIMENTALLY OB- 

SERVED VALUES (IN PARENTHESES) 
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Actinide compounds have Important roles In nuclear fuels, and 

recent attempts to develop more efficient methods to enrich uranium using 

lasers have focused attention on UFg. In this connection, we carried out 

calculations on UFg to expand on the fragmentary knowledge about the nature 

of chemical bonding In actinide compounds. We Illustrate here how some 

of the calculated properties compared with their experimentally observed 

values. The bond length of UFg In the ground state configuration was 

( 31 1 

optimized. The predicted value of 3.70 bohr was found In good agreement^ * 

with the experimental value of 3.777 Bohr. 

In Table 7 we compare the calculated excitation energies In UFg', 

after Introducing sp1n>orb1t coupling using the simple approach defined by 

Equation (31), with the experimentally determined energies. The comparison 

is quite favorable, with the theoretical values only 0.1-0.2 eV higher In 
(31 1 

energy.' ^ Table 8 compares experimental vertical Ionization potentials of 
UFg with orbital energies calculated using Koopman's theorem. The orbital 
energies shown In this table overestimate the experimental Ionization poten- 
tials by nearly 3.9 eV. Shifting the calculated values by 3.9 eV brings the 
calculated levels Into harmony with the experimental photoelectron spectrum. 
Similar experiences have been found In all electron ab Initio calculations 
on SFg, for example. The problem Is well known; it lies In the Inadequacy 
of the delocalized molecular orbitals of the neutral molecule for describing 
the relaxation effects associated with the "localized hole states" of the 
Ionized molecule. Since relaxation effects appear to alter only the 
absolute positions of the states of UFg^, an assignment could be made of the 


iiotoelectron spectrum by use of the shifted results In Table 8. The first 


Table 7. Excited States of UFg . Comparison of Theoretical and Experimental Results 


(Ref. 31) 




Excitation energy (eV) 


State 

Exptl 
(Bef. • > 

Bel. ECP 

(R«f. 31 ) 

DVM 
(Bef. b ) 

Xa-SO 
(Bef. c ) 

r,w 

0.00 

0.00 

0.00 

0.00 

Ta, 

0.57 

0.67 

1.00 

0.85 

r». 

0.86 

0.97 

1.13 

1.06 

Ta, 

1.58 

1.80 

2.51 

2.53 

I'a* 

1.77 

1.95 

2.69 

2.50 


J. Reisfeld and G. A. Crosby. Inorg. Chem. • 4, (1965). 

D. Koelling. D. E. Ellis, and R. J. Bartlett, J. Chem. Phys., 65, 3331 (1976). 
^A. N. Boring and J. H. Wood. 
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Table 8. Comparison of Experimental Vertical Ionization Potentials of UF- 

D 

with calculated orbital energies (Koopmans' theorem), (ref. 31) 



Bel. HF 

Rel. HF with 8pin-orbit coupling 

Exptl* 

«1. 

18.43 

8k 

17.88 

(13. 98)* 

13.9 


Su 

19.11 

(15.21) 

15.35 


18.56 


18.55 

(14.65) 

14.8 



Sg 

18.57 

(14.67) 



19.64 

lu 

19.63 

(15.73) 

15.98 


8k 

19.65 

(15.75) 


3a 

19.73 


19.73 

(15. 83) 



20.11 

8k 

20.09 

(16. 19) 

16.58 


6k 

20.16 

(16. 26) 


ih, 

20.58 


20.55 

(16.65) 

16.85 




20.59 

(16. 69) 


2^/ 

20.74 

^g 

20.74 

(16. 84) 

17.30 

2fj„ 

37.76 

8k 

35.20 




6k 

41.95 



2® If 

44.50 

^g 

44.50 



U, 

44.71 


44.71 



Iti, 

46.44 

8k 

46.04 




6k 

48. 56 



Uu 

70. 92 


70. 92 



Jrel. hf with spin-orbit coupling results shifted by -3.9 ev. 

KARLSSON, L. MATTSSON, R. JODRNY, T. BERGMARK. AND K. SIEGBAHN. 
PHYS. SCR. 14. 230 (1976). 
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peaks then coincide, and the remaining peaks lie no farther than 0.6 eV 
from a predicted Ionization potential. 

5. JJYNOPSIS 


We close this part of this paper by listing the main character- 
istics of the method discussed here. 

(1) This approach to the calculation of molecular relativistic 
electronic structure Is based on one-component wavefunc- 
tlons familiar from the nonrelativistic theory. The two 
principal Ingredients to this approach are: (a) Cowan- 
Griffin relativistic Hartree-Fock atomic orbltals^®^, 
and, (b) "relativistic" effective core potentials. 

(2) The relativistic effects on molecular properties of 
chemical Interest such as the relativistic effects on the 
bond lengths and on the energy ordering of states appear 
from all cases examined to be reliably obtained by this 
method. The molecular wavefunctlons and energies are first 
calculated Including only the relativistic mass-velocity 
and Darwin terms. If the effect of spin-orbit coupling 

Is expected to be Important, then the results of the first 
calculation are combined with the calculation of the spin- 
orbit matrix. The latter Is calculated using an effective 
one-electron spin-orbit Interaction approximation. The net 
Interaction matrix obtained In the end Is diagonalized to 
obtain the final molecular states. 
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(3) This approach allows the molecular cal itlons to be 
perfoimed by direct application of the tull range of 
the traditional quantum chemistry techniques. Including 
NCSCF and Cl, to Include electron correlation In ob- 
taining the ground and excited states of molecules. All 
computational aspects of the approach that pertain to the 
RECP's have been tested and Implemented In standard quan- 
tum chemistry programs. 

(4) This approach Is applicable to polyatomic molecules of 
arbitrary geometry. Gradient assisted RECP searches 
for the critical points of the potential energy surfaces of 
a transition metal complex have recently been reported. 

(5) Large computational simplifications over traditional all- 
electron treatments are obtained via the "relativistic” 
effective core potentials In their role as a device to 
reduce the molecular problem to just the valence electrons. 

6. CRITICAL ASPECTS OF THE CONSTRUCTION OF RECP's 

The reliability of the RECP's for molecular calculations rests 
on the properties given to the valence pseudo-orbital through Its 
definition, and the details thereof are thus of critical Importance. Indeed, 
every Improvement In the approach has been associated with a revision of the 
properties of the valence pseudo-orbital. 
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5.1. Linear Combination of RHF Orbitals 

The traditional definition of a valence pseudo-orbltal has been 
to take It first as a linear combination of the core and valence orbitals 
of the same orbital angular momentum, 

bttl 

The coefficients were chosen to satisfy the broad criteria given in 

Section 2.2. These criteria were indeed satisfied almost optimally by 
requiring, in the following order, that:^^^^ 



(b) P te normalized 

4 

M 

(c) P maximize a functional of .the "smoothness" of the orbital. 

‘al 

The major deficiency of this definition has been found to lie in the restric- 
tion of the pseudo-orbital to the expansion given by Equation (32). Since 

% 

the valence pseudo-orbital is normalized, the expansion by necessity has 
less amplitude than the original valence orbital at large radial distances. 
The net effect is to introduce a charge redistribution into the valence space 
that is unphysical and can seriously impair the reliability of calculated 
potential energy curves. A comparison between the Hartree-Fock valence 
orbital and the valence pseudo-orbital is shown in Figure 14 that illustrates 
the charge redistribution for the case of the silicon atom 3s orbital. 


AMPLITUDE 
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Silicon "I 
3s orbital 


Li near-Combi nati on-of- 

Atomlc-Orbitals 

Pseudo-orbital 


"Shape-Consistent" Pseudo-orbital 

/ 



Radial Distance, Bohr 

FIGURE 14. COMPARISON OF THE HARTREE-FOCK VALENCE ORBITALS WITH THE VALENCE 
PSEUDO-ORBITALS OF THE "LINEAR-COMBINATION-OF-ATOMIC-ORBITALS" 
APPROACH AND THE "SHAPE CONSISTENT^' APPROACH FOR THE CASE OF THE 
3s ORBITAL OF SILICON ATOM IN THE ’’P STATE 
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6.2. Shape-Consistent RHF Pseudo-Orbitals 

The remedly to this deficiency has been found to lie In abandoning 
the traditional expansion* Equation (32)* as a part of the valence pseudo- 
orbital definition. Instead* motivated by the Importance of preserving the 
charge distribution characteristics of the original valence RHF space* the 
valence pseudo-orbital Is first defined as Identical to the WF valence 
orbital over the valence segment of the orbital 


K 





(33) 


where 1s a radial distance larger than the outermost node of P^^* yet 
not much larger than the radius of the outermost valence maximum. The defl- 
nitlon of In the core segment of the orbital* 0 ^ r ^ Is based In 
turn on satisfying at least the requirements of: 

(a) Nodelessness 

(b) Matching of end-point (at r*0 and r»R|^) continuity 
conditions on the orbital and Its derivatives (up 
to the third derivative at r=Rj^) 

(c) Normalization. 

The core segment of P^^j^ Is expanded In polynomials^ '* or even exponential 
functions.^ However, only as many parameters are Introduced In these 
expansions In practice as are required to fulfill the above listed require- 
ments. The criterion of smoothness has not been used explicitly In this 
approach. Only the choice of R|^ has loosely been connected In practice 
with the desired ultimate smoothness of P^j^*' ' 
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The valence pseudo>orb1ta1s obtained by this definition are 
referred to as “shape-consistent" pseudo-orbitals. A typical "shape- 
consistent" valence pseudo-orbital Is compared with the previous pseudo- 
orbital obtained from the llnear-comblnatlon-of-atomic-orbitals defini- 
tion In Figure 14 for the case of the silicon atom 3s orbital. The 
restoration of the proper amplitude In the valence region Is shown In 
Figure 14 to be accwnpanled by a comnensurate adjustment In the amplitude 
over the core region. Whereas the "shape-consistent" definition of valence 
pseudo-orbitals has been discussed by a number of workers, Christiansen, 

Lee, and Pitzer were the first to demonstrate the Implications of the "shape- 
consistent" ECP's for the reliable calculation of potential energy curves. 

One of their conclusive results for the case of the potential energy curve 
of the CI 2 X'lg^ molecule Is shown In Figure 15. The result of the "shape- 
consistent" ECP Is shown In this figure to be In excellent agreement with 
the comparable all-electron (AE) result. In contrast, the previous ECP 
clearly led to a much poorer comparison. Figure 15 also shows another 
potential energy curve In reasonable agreement with the all -electron result. 

This other potential energy curve was obtained by an Intermediate remedy 
(42l 

to the problem' ‘ that has now been superseded by "shape-consistent" 
approach. 


6.3. Hamiltonian and Shape Consistent RHF Pseudo-Orbitals 

The replacement of the RHF valence orbitals by the valence pseudo- 
orbitals In the valence- valence Interactions causes a modification In the 
Interactions commensurate to the modification In the orbitals. There Is a 
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R (bohr) 

FIGURE 15. COMPARISON OF THE ALL-ELECTRON POTENTIAL ENERGY CURVE FOR THE 
Cl? MOLECULE IN THE X'ri STATE (CALCULATED WITH A GENERALIZED 
VALENCE BOND WAVEFUNCTIDN) WITH THE CORRESPONDING VALENCE 
ELECTRON RESULTS USING ECP's FROM THREE DIFFERENT APPROACHES: 
THE "SHAPE-CONSISTENT" APPROACH (CHRISTIANSEN, LEE, AND PITZER; 
REFERENCE 39), THE OLD "LINEAR-COMBINATION-OF-ATOMIC-ORBITALS" 
APPROACH. AND THE "INTERMEDIATE-FIX" APPROACH. 
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change In these Important Interactions even In the case of the "shape I 

consistent" definition of the valence pseudo-orbitals. In this case, 1 

J 

although the "shape-consistent" psuedo-orbital Is Identical over the valence 

segment with the RHF orbital, the core segment of the pseudo-orbital only j 

satisfies the normalization condition among all the Important other orbital 

properties that determine the valence- valence Interactions. Rappe, Smedley, 

and Goddard have proposed to enhance the definition of the core segment of 

the valence pseudo-orbital so as, to minimize the remaining error In the 

( 43 ) 

valence- valence Interactions.' ^ This enhanced definition of the valence 
pseudo-orbitals Is referred to as the "hamlltonlan and shape consistent" 
approach. 

6.4. The Method of Moment- Accumulation Functions: A New Approach 

The minimization of the error In the valence-valence Interactions 
supplies Important additional conditions for the definition of the core 
segment of the valence pseudo-orbital. However, the relationship between 
the satisfaction of this criterion and the requisite characteristics of 
the Individual core segments of the valence pseudo-orbitals Is complex and 
Indirect. We show below a new approach whereby the specific orbital proper- 
ties that affect the minimization of the error In the valence-valence Inter- 
actions can be Isolated from the peripheral Intricacies of the energy Inter- 
actions.' ^ In fact, these new specific orbital properties are found to 
constitute a set of conditions on the Individual orbitals of the same 

(44) 

simplicity as the normalization condition.' ’ 



r 


The sum of the valence- valence Interactions Is contained In the 

(45) 

atonic valence energy. The atomic valence energy Is given by' 
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The terms In this expression range only over the set of valence orbitals 

denoted as )[. The orbital occupation numbers are denoted while 
the A (£,£') coefficients are the Roothaan vector coupling coefflclents^^®^; 

A 

the latter are simply related to the Clebsh-Gordan coefficients. The 
terms are one-electron Integrals of the effective one-electron operator 
containing the RECP. The term stands for the deviation of the multi pi et 

energy from the average energy of configurations. This energy contribution 
Introduces only additional combinations of the type of Interactions already 
present, and Its detailed structure Is therefore of no Importance here. 

The two-electrons Interactions among the atomic valence orbitals 


are given by the well known F and G integrals. These Interactions are 
,(45) 


defined by' 


F V*T) « I f F ,W it it' 


and 
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In this traditional form, the properties of the core segment of the valence 
pseudo-orbital are Interwoven In the Intricacies of the energy expression, 

E'”'', and no discrete succinct measures of the valence pseudo-orbital 
characteristic can be seen to relate to the minimization of the error In 
the valence- valence Interactions. 

We show next a new expression for the formula of the F and G 
Integrals that does Indeed bear out through these Integrals the dependence of 

the valence-valence Interaction on a set of discrete succinct measures of orbital 
characteristics. We define the function 

i r I , 

We shall name this orbital characteristic the “moment accumulation" function. 
These functions become the various x-order transition moments of the valence 
set of functions as the radial distance from the nucleus approaches Infinity. 

The special case of the function describes the physical Increase In 

the amount of charge deriving from the orbital as the radial distance 

from the nucleus Increases; we refer to It as "charge accumulation" function. 

In particular 


twv. Q* (<•) = I 


(37) 



because of the normalization of the radial functions. Analogously. 

•lUw Q tr) - 0 (38) 

f^co wl,WI 

because of the orthonormallty of the radial functions. 

The function Q® _,(r) first serves to clearly bring out the 
nature of the Improvement of the orbital properties brought about by the 
'thape consistent" pseudo-orbital definition. Figure 16 shows a comparison 
of the Q 3 j^ 3 s(»*) function deriving from the Hartree-Fock valence orbitals 
with those deriving from the valence pseudo-orbitals for the case of the 3s 
orbital of silicon atom. Apart from an expected small discrepancy In the 
core region, the "shape consistent" 3s pseudo-orbital leads to an exact match 
with the Hartree-Fock charge distribution over the whole valence region. The 
overall Improvement In the charge distribution over that obtained from the 
older pseudo-orbital deriving from the llnear-comblnatlon-of-atomic-orbitals 
definition Is shown In Figure 16. The latter pseudo-orbital yields a charge 
distribution that differs from the Hartree-Fock distribution not only over 
the core region but over a large part of the valence region. 

We note, as a digression, that the "charge accumulation" function, 
3 ^(r), serves also to give a measure of the onset of the physically meaning- 
ful valence segment of a valence orbital. The outermost maximum of the 3s 
orbital Is shown In Figure 14 to be at about r«1.26 Bohr. This point might 
perhaps be Interpreted as already In the midst of what one would call the 


CHARGE ACCUMULATION 
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FIGURE 16. COMPARISON OF THE CHARACTERISTICS OF THE "CHARGE ACCUMULATION" 
FUNCTION OF THE HARTREE-FOCK 3s ORBITAL OF SILICON ATOM, AS A 
FUNCTION OF RADIAL DISTANCE, WITH THE "CHARGE ACCUMULATION" 
FUNCTION OBTAINED FROM THE VALENCE PSEUDO-ORBITALS OF THE "SHAPE- 
CONSISTENT" APPROACH AND THE OLD "LINEAR-COMBINATION-OF-ATOMIC- 
ORBITALS" APPROACH 
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physically meaningful valence region. We note, however, that Figure 16 
shows the amount of charge accumulated up to this outermost maximum, at 
r>1.26 Bohr, to be only about .05, 5t of the total charge. 

The properties of the “charge accumulation" function are also 
directly related to the basic electrostatic Interactions among the electrons. 

The traditional expression for the Coulomb potential with the 3s orbital density 
as Its source Is 



Carrying out the angular Integrations yields the well known formula 



(40) 


The contribution of the orbital Is spread out over two distinct terms 
In this expression. Integrating by parts, however, we find that both terms 
can be consolidated Into a single term that depends concisely on 
through the "charge accumulation" function, «Ss.3s<-). We obtain 
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The asymptotic expansion of the above expression for the Coulomb potential 
for large r yields the well-known poln^ charge-like expression 
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( 42 ) 


This Illustrates the central role that the *rb1ta1 properties contained In 
the charge accumulation function have in the basic electrostatic interactions 
among electrons. 

Figure 17 shows a comparison of the Coulomb potentials resulting 
from the three different "charge accumulation" functions shown previously In 
Figure 16. The Coulomb potential deriving from the "shape-consistent" 
valence pseudo-orbital Is seen to agree perfectly In the valence region with 
the potential derived from the Hartree-Fock orbital. In contrast, the 
pseudo-orbital deriving from the 11near-comb1nat1on-of-atom1c-orb1ta1s 
definition yields a Coulomb potential that, consistent with the corresponding 
"charge-accumulation" function shown In Figure 16, Is too repulsive even In 
the valence region of the valence orbital. 

Indeed, the F and G Integrals themselves may be reexpressed In 
terms of the "moment accumulation" functions. Starting from the definition 
given by Equations (35a) and (35b), and Integrating by parts repeatedly, 
we find that 


I 

fVwi') = ltV«o] 


(«•) 


and 
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FIGURE 17. 


COMPARISON OF THE CHARACTERISTICS OF THE COULOMB POTENTIAL. AS 
A FUNCTION OF RADIAL DISTANCE. OBTAINED FROM THE DENSITY OF THE 
HARTREE-FOCK 3s ORBITAL OF SILICON ATOM WITH THE COULOMB POTEN- 
TIAL OBTAINED FROM THE VALENCE PSEUDO-ORBITALS OF THE "SHAPE 
CONSISTENT" APPROACH AND THE OLD "LINEAR-COMBINTATION-OF-ATOMIC- 




56 OF 

This proves that through these expressions for the F and 6 Integrals* the 
valence- valence Interaction In E . Equation (34)* depends on the valence 
pseudo-orbital properties embodied in the set of "moment accumulation" 
functions. 

It follows that minimization of the error In the valence- valence 
Interaction Is equivalent to the requirement that the core-segment of the 
valence pseudo-orbitals In the "shape consistent" definition satisfy the 
Hartree-Fock values of 

and 

I (fc,) , I =«'*’>. » I . i 

^ ^ (44b) 




This set of conditions contains the normalization conditions* 

as a subset. The replacement of the single condition of normalization by 
the above "moment accumulation" set of conditions appears to provide the 
most natural transition from the "shape-consistent" Into the "hamlltonlan 
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and shape consistent" approach. Indeed, the increase of the required set 
of "moment accumulation" functions beyond Equation (44) provides a systematic 
procedure ■*^or the definition of pseudo-orbitals that increasingly tend to the 
original Hartree-Fock valence orbitals and is independent of the availability 
of a pertinent energy expression. The limit is, of course, inconsistent with 
the desired criterirn of nodelessness in the definitions of valence pseudo- 
orbitals, and is, therefore, primarily of interest as a generalization of the 
approach. 

The transition from the "normalization" condition to the "moment 
accumulation" conditions in the "shape-consistent" approach requires an 
extension of the form of the expansion of the core segment of the valence 
pseudo-orbital. Rather than simply expand in a single continuous polynomial 
of increasingly higher degree, it appears to be preferable to expand in a 
large number (p''ssibly simulating a complete basis) of piece-wise continuous 
polynomials such as, for example, cubic splines of finite support. Moreover, 
the condition of "smoothness" of the valence pseudo-orbital can simply be 
added to determine all degrees of freedom in the expansion not fixed by the 
"moment-accumulation" conditions. This provides a systematic procedure 
for obtaining valence pseudo-orbitals wherein the final form of the resulting 
RECP's is independent of constraints on its shape that may result from an 
expression fo»^ the core segment of the valance pseudo-orbital with just 
a minimal set of parameters. Finally, we note that the "moment accumulation" 
conditions include conditions that couple the various valence pseudo-orbitals. 
The satisfaction of the "moment accumulation" conditions hence requires that 
all the valence pseudo-orbitals be solved simultaneously. In contrast, the 
old "shape-consistent" procedure allowed each valence pseudo-orbital to be 
obtained independently since no recognition was given to the interdependence 
of these orbitals. 
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7. THE L0N6-RAN6E BEHAVIOR OF THE RECP: A NEW ANALYSIS 


The characteristics of the long-range segment of the RECP's have 
been shown to be an Important factor In the reliability with which molecular 
potential energy curves can be calculated using the RECP's. This Inter- 
dependence Is not surprising In view of the large overlap between the long- 
range segment o^ the RECP's and the electron densities of neighboring atoms 
at equilibrium Internuclear distances. Atomic calculation: , In contrast, 
provide almost no test of the long-range properties of the RECP's because 
of the small overlap of the atomic charge density with the long-range se^nt 
of the RECP. 

We present next the first detailed analysis In terms of moment 
accumulation functions of the characteristics of the long-range behavior 
of RECP's. Starting from the expression for the RECP given In Equation (18), 
and using the properties of the P^^^^ orbitals, one obtains 



for r ^ Rj^. The first term Is the local core-valence Interaction, and the 

second tenn Is the local residual In the valence- valence Interactions. These 

terms consist of radial Coulomb and exchange potentials, the general form 

of which Is traditionally given In terms of the Y functions. The latter are 
defined as 
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We first obtain an alternative expression for the Y functions in terms of 
“moment accumulation" functions. 
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(48) 


by Integration by parts. 

We first analyze the residual of the valence* valence Interactions. 
Using the expressions for the Coulomb and exchange potentials in terms of 
the "moment accumulation" function, we obtain 


ifoi 




Vjfev 




(49) 


The residual of the valence- valence interactions Is seen to depend on the 
difference between the "nxxnent accumulation" functions calculated from the 
valence RHF orbitals and from the valence pseudo-orbitals (indicated by a 
superscript tilde). We have restricted the expression for the valence- valence 
interaction to the case of the average-energy-of-configurations since it 
suffices to illustrate the long-range characteristics of the RECP. 


ORIGINAL Fk>: - *3 
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The core segment of the valence pseudo-orbital of the "shape 
consistent" approach satisfies the normalization condition 


Consequently, In the “shape-consistent" approach the first and third terms 
of the residual In the valence- valence Interaction, Equation (49), are zero. 
However, the remaining terms In the residual of the valence- valence Inter- 
actions are not equal to zero for “shape cons1stent"-only orbitals. In the 
"moment accumulation" formulation of the "hamlltonlan and shape consistent" 
approach, we require that the core segment of the valence pseudo-orbital 
satisfy the full set of conditions 




Y.!eV> W/'«y 

(51b) 


It follows from Equation (49) that in the "hamiltonian and shape consistent" 
approach the full residual in the valence-valence interactions vanishes for 
r ^ Rjy|, i.e. , 
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The RECP is therefore reduced to the local core-valence interaction. 


U^CD 
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(53) 


provided the moment accumulation conditions given by Equation (51) are satis- 
fied. Expressing the Coulomb and exchange interactions in terms of "moment 
accumulation" functions, we obtain for r ^ 



/ <D i ^ 1 


(54) 


The set of all core orbitals is denoted as Jj^. 

Finally, we discuss the long-range behavior for r»R|y|. The 
asymptotic expansion of the Y functions for large radial distances yields 
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In turn* the asymptotic expansion of the RECP for large radial distances Is 


itf 






P.e i 


(56) 


under the assumption that the residual In the valence- valence Interactions 
Is zero. The leading term Is the point-charge potential Nj./r where Is 
the number of core electrons. This Is Indeed In accord with what Is expected 
on physical grounds. If the core segment of the valence pseudo-orbital only 
satisfies the normalization condition, then the residual In the valence- 
valence Interactions unfortunately also contributes a term to the long-range 
Interaction. We find In that case that, for r»Rn^, 



This long-range contribution, however. Is clearly an artifact of the 
Insufficient definition of the valence pseudo-orbital. 
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Figure 8. 


Comparison of the non-relativistic (NR) and relativistic (R) 
Hartree-Fock 7s orbitals In uranium atom. 

Comparison of the non-relativistic (NR) and relativistic (R) 
Hartree*Fock 5f orbitals In uranium atom. 

Pi 

A nodeless 2s orbital In lithium atom, , and the corre- 

2s 

spending nodeful 2s Hartree-Fock orbital. 

The s, p, and d local effective core potentials for the 
lithium atom. 

Schematic representation of the one-electron orbital spectrum 
associated with the s, p. and d local effective core potentials 
for the Lithium atom, respectively. 

Comparison of a valence molecular orbital for LiH X'r‘*‘ obtained 
from the valence electron (VE) ECP calculation and a comparable 
all-electron (AE) calculation. 

4. 2 + 

Comparison of the valence molecular orbital for Li 2 X 
obtained from the valence electron (VE) ECP calculation and a 
comparable all-electron (AE) calculation. 

% 

Schematic comparison of two possible valence pseudo-orbitals, P^^^^ , 
with the parent Hartree-Fock valence orbital, P^^^. 

Schematic comparison of two possible valence pseudo-orbitals, P^^. 


Figure 9. Comparison of the nonrelativistic (NR) and relativistic (R) 

valence electron (VE) potential energy curves for the low-lying 
electronic states cf the XeF molecule prior to the inclusion 
of the spin-orbit effects. 
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Figure 10. 


Figure 11. 


Figure 12. 


Figure 13. 


Figure 14. 


Figure 15. 


The relativistic (R) valence electron (VE) potential energy 
curves for the low>ly1ng electronic states of the XeF molecule 
Including spin-orbit effects (using the atoms- In-molecules 
approximation). 

Spin-orbit coupling matrix elements between electronic states 
of the rare gas oxides computed by Langhoff (Reference 22) 
using the full Brelt-Paull hamlltonlan (solid lines) and the 1- 
electron 1 -center effective spin-orbit operator (circles and 
crosses) 

The potential energy curves for the ground states of AuH and 
AuCl obtained using a nonrelativistic (NR) and relativistic (R) 
ECP and a generalized valence bond wavefunctlon. 

Geometrical and energetic parameters calculated by Hay (Reference 
30) for Zelse's salt, and comparison with experimentally ob- 
served values (In parentheses). 

Comparison of the Hartree-Fock valence orbitals with the valence 
pseudo-orbitals of the "llnear-comblnatlon-of-atomic-orbitals" 
approach and the "shape consistent" approach for the case of 

3 

the 3s orbital of silicon atom In the P state. 

Comparison of the all -electron potential energy curve for the 
Cl 2 molecule In the X'l^ state (calculated with a generalized 
valence bond wavefunctlon) with the corresponding valence 
electron results using ECP's from three different approaches: 
the "shape-consistent" approach (Christiansen, Lee, and Pitzer; 
Reference 39), the old "llnear-comblnatlon-of-atomic-orbitals" 
approach, and the "Intermediate-fix" approach. 
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Figure 16. Comparison of the characteristics of the "charge accumulation" 
function of the Hartree-Fock 3s orbital of silicon atom, as a 
function of radial distance, with the "charge accumulation" 
function obtained from the valence pseudo-orbitals of the 
"shape-consistent" approach and the old "linear-combination- 
of-atomic-orbitals" approach. 

Figure 17. Comparison of the characteristics of the Coulomb potential, as 
a function of radial distance, obtained from the density of the 
Hartree-Fock 3s orbital of silicon atom with the Coulomb poten- 
tial obtained from the valence pseudo-orbitals of the "shape 
consistent" approach and the old "linear-combination-of-atomic- 
orbital" approach. 
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PROGRAM SHELLIINRUT ; 3UTPUTI 

c ■ ■ ■ 
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?10 CONTINUE 
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t.THlS ROUriNt EVALUATES THE ANGuLM HONENTUN COUPLING CONSTANTS 
;OH«Or;/FTABCH/FPQRU3t t3fl5l 

CONHON/ZTABC M/LKF U3) »LNL (491 • ZLH (130) »LHX (1901 tLMVT 130 tLNZ (130 > 
IOsL*(L*l)-N+l 

IMN«LNF(IC) 

1H*»LNL(10) ‘ 

J0*LP*(LP*1)-NP*1 
JHN«LNF( JD) 

JNXsLNLCJOl 

■ SUMIsO.OEO ~ 

CO ZC «:1HN,1MX 

SUMJsO.OEO 

00 IQ JsJMM.JHX 

JX = LMX (I) ♦KX+LMX (J) »1 

_ J]^LNy III*KY*LHJ[( Jl+l 

JZ*LHZ(ii*KZ»LMZ(j)*l 
10 SOMJsSUNJ ♦ Zi.r: J)«rpQRUXt JYt JZ> 
zr SUMIsSUHI ♦ ZtN(II*5JNJ 


ORic:^;>tu t-. . - ;s 
OF POOR QUALITY 


FUNCTION tCO<KX,KYt<ZtLtH tKX P,XyPtKZP> 

CONNON/7TABCH/LHFU9) flNLU?) tZLNIlSO) *LHX 1130 1 «LMV( 1301 1 LHZ 1 1 ^t0» 
UOHNON/FTABUH/FPQRdS. 13il3) 

lost* a»ll-M»l 

IHN«i.HFaO> _ _ 

INX>i.HL(10) 

SUH«u.Q£0 
00 10 Z>INN,IHX 
AXsXX*LNXtU«<XP*l 
lY»<»*LHYtll AKYPAl 

iZsKZ»LHZiii*kZP*i ' " 

SUNsSUHAZLNIil •FPQRUXilYtlZ) 

10 CONTINUE 

EC0«SUN 

RETJ^N 

ENO 


OF POC- ^ hY 


SUBROUTINE FTIB _ 

C...TH1 S Routine SEts up i table of F*FUN;ffON Values 

c 

INTEGER PfQtR 

C0HH0N/FTA6CM/F>QR(139 13tlSI 
OATA PI/ 3.1<ii592EB3B69SE0/ 

C 

C*,. ZERO-OJT the fABLE 


DO 30 P*lfl3 
DO 2D Q«l«13 
DO 10 R«itl3 
10 FPQ2(P,atP>>0«0E0 

20 CONTINUE 

30 CONTINUE 

C 

C... RECURSIVELY GENERATE THE NON-ZERO ENTRIES TO THE TABLE 
FPQR«l,l,l)*ii. 0E0»PI 
DO aO P*lfl3#2 
r P*P-1 

DO 70 Q=iti3,2 

QQsQ-1 

00 60 R«ltl3,2 
RR»R-1 


C 


IF IP.EO.ll GO TO 

FPQRCP*Q,RI* (PP-l, dEai ♦FFQR'(P-2,Q,RI /I PP«^QG^RR • CEO I 
GO TO 60 

40 IF (Q.EQ.ll GO TO 90 

FPQRIPtQtRMfQQ-l.OEOMFPQRIP, Q-2» R ) / C PP ♦QQ»RR4l. 6E0 ) 
GO TO 60 

SO IF IR.EQ.l! GO TO 60 

FM? IPt(j,R»»tRR-i.CE8>*F PQR (P* Q,e*2) / IPP»QQ»RR«1.0E0 ) 


60 CONTINUE 
7C CONTINUE 
SO CONTINUE 


C 


ketU%N 

END 




C*''' »Li'Ty 


SUBROUTINE ZT»B 

C...TH1S routine sets up fNE'REAL SPHERlClL HRRHONiSS'iM The ToSm“ 0? LINERR 
(....COMBINATIONS OF CARTESIAN PRODUCTS-lLtNIsLILHi -N»l 

COHNON/Zr ABCN/LNFU9) tLNLI<i9l»ZLN(lS3l*LHX(t30) tLNrilSOI .LHZI130I 

C 

DATA FFIF12.S663r06t%3$9E0/ 

DATA LHF/l«2t3.%t9tZ.Btl0t il.lZf 19*16»i6*2^22«23»Z5.2^3a.3A.S6t 

X39. Al* A3, 4$, (.7, 90,53, 92/61. '6%* 67 ,707r2,Z6«T». 61* «<r. 07.93,97.10 3. ' 
X1Q6, 110, 11 3, 11 6, 120, 124,127/ 

CATA CM./1, 2, 3, A, 6, 7, 9, 10, 11, 13, 19, 17. 19, 21, 22, 24,27, 29, 33,39,30, 

XAO, 42, 44, 46, 49 ,92, 56, 60, 63, 66, 69, 71, 75 ,77, SO ,03,06,92,96,102, 1C9, 
X109, 112, 119, 119,123, 126, 130/ 

DATA I.NX/0,l,0,0,2,a,l.,0,0,0,l,3,i,2,0,l,l,0,0,0.0,l,2,0,4^2tO,3, 

XI, 2, 0,0, 2, 1,1, 0,0, 0,0, 0,1, 1,2,6, 3, 1, 5, 3,1,6, 2,4,3,1, }.,3,2,0,2, 6, i, 

XI, 1,0, 0,0, 0,0, 0,1, 1,2, 2, 0,0, 3, 1,4, 2, 0,9, 3, 1,5, 3, 1,4, 2, 0,4, 2, 0,3,1, 

XI, 3, 2, 4, 2, 0,2, 0,1, 1,1, 0,0, 3, 0,0, 0,0, 1,1, 1,2, 3, 2, 0,3, 1,3, 1,4, 2, 0,6, 
X4,2,0/ 

DATA LNY/0, 3, 0,1, 0,2, 0,0, 0,1, 1,0. 2, 0,2, 0,0, 0,0, 1,1, 1,1,3, 0,2, 4,0, 

X2, 0, 2, 2, C,C,0, 0,0, 0,1, 1,^,1, 1,3, 1,3, 0,2, *,,4, 2, 0,0, 2, 2, 0,0, 2, 0,2,0, 
X0,6,0,G,U,i,l,l,i,l,i,l,3,3,l,~3,l,3,9,l,3,9,0,2i4,0,2^’',9,2,4,0,2, 

X2, 0,0, 2, 0,2, 0,2, 0,0, 0,0, 0,0, 0,1, 1,1, 1,1, 1,1, 3, 1,3, 1,3,1, 3, 1,3, 5,0, 
X2,4,6/ 

DATA L HZ/(, , 0, 1, 0, 0, 0, 1, 2, 0,1,0,(^, 0,1, 1,2, 0,3, 1,2, 0,1, 0,0,0, 0,0,1, 

XI, 2, 2, 0,0, 3, 1,4, 2, 0,3, 1,2, 0,1,1, 0,0, 0*0, 0,1, 1,1, 2, 2, 0,0, 3, 3, 1,1, 4, 
X2,C,^, 3,1,4,2,4, 3, 1,2,0 ,_2^0 1 ,0,0 , ^ 0,0, 1,1, 1,2, 2, 2,0 ,0,0, 3,3, 
Xl*l,4,4,2,2,0,6*9,3,i,E*4,2,0,$, 3, 1,4,2,C*3, 3,i,l,2, 2, 1,1, 1,0* 
Xb,0,0/ 

C 

ZLH(l)ssaRT(1.0E0/FPn 

ZLH(2lx3gRT(3.0E0/FPl) 

ZLMIJI*ZtHI2) 

ZLMC4>»ZLH(2r - -- - 

ZLM(9)xSQRT(19.0E3/FPn/2.0EO 

ZLH(6lx-ZLMI9) 

ZLH»7)*2.0EI)*Zl.H«5» 

ZLN(Blx3.0E0*SQRT (9. dEO/FPl ) /2 .DEO 

ZLH(9Ix-ZLHI8>/3.0E0 

ZLNIlOlsZLHiri 

ZLH(11>xZlMI7> 

ZLH(12)xSQRT(39.OEO/(6.OE0*FPnt 

ZlMll3)x-3.3Ei«*ZLHtl2) 

ZI.H(14IxSORT( 139.0C0^(4. 0E0*FPIt I 
ZLHU5l»-ZLHa4» 

ZLMU6»»5.0E0*SaRTi2i;eE0/tS',0E0*FPrM 
ZtHll7)x.Zi.H(16)/9.0£0 
ZLMa8l*5.0£0*SQR^ir,0E0/FPI»/2. OiO 
ZLN(19l«-3.aE3*ZLH(19»/9.0E0 
ZLH(20)xZLMI 16) 

2LM«2l>x2LM(ir» 

ZLHI22Is2.6E0*ZlHI14» 

ZLrt(23»«*ZLM(n» 

ZI.>M24I«-ZLH(12I 

ZLH(29)«SQRT(319.CE0/<64.0EC*FPII» 

ZLHI26>*-6.0E0*ZuHI25» 

ZLM(27|x7umi25j 

ZLMI28»xSQivT«315.0E3/ 18. 0E0»FPin 
ZLN(29lx.J,9tO>ZLHI29) 


I 


*■ ■ 

OF 


TEHP>-SQRT US.OE0/FPI)/4^OEO 

ZkHI36»«7.CE0*nHP 

ZLH(3t>»ZLHI30) 

ZLHI3Z)>TEHP 

ZLni33)«-TEHP 

TEHP«SQ?T U5.0E0/lt.0EO«FPlM 

ZLM34>aZ.0EC*TEHP 

ZLKI3S)a>3.6E6»fEHP 

TEHPsSOPT(9.0EOZFPl>/a.OEO 

ZLHI36»35.0E0*TEHP 

ZLMI37)»-3u.0ED*TEM» 

ZLHt38»>3.0E0*TEMP 

ZLHI39)xZLHI34> 

Zi.H(40i*ZLH(39) 

TENP«SORM45.0E0/U.eE0*FPIM 

ZLM|I»1J«Z.0E0*TEHP 

ZLHUZta'TEKP 

ZLHI43)s-ZLHtZ9) 

ZLH|44)«-ZLHIZ8I 

ZLHi45|sSQKT(3i5.0£3r|4. OEO*FPI»» 
ZLHI46I>-ZLHI45I 

ZLHU7}sSQRTI693.0E0/(129.0E0aFPI) I 
ZLHI4S fa-10. CEO* ZLN( 47 > 
ZLm49ls5.0E0*ZLHI47f 
ZLHtSOf«SQRT(3469.(l£0/l64.0EO«FPl) I 
ZLM(91)*-6.9E0*ZlN(90I 
ZLN($2}«ZLH(S0t 

TEMP>SORT(39S.aiO/ll20.ilE9*FPIfl 
ZLN(53)s9.0EQ*TEHP 
ZLMI54)«-27.0E0*TEMP 
ZLMI551«3.0EC*T£HP 
ZLMj50)»-tEMP ■ 

TEMPiSORT (115$.SE07FPI>/4.0EO 

ZLHf57»a3.6£0*T£HP 

ZLH(S0)*-ZIH(S7) 

ZLM|59>»-TE*1P 

ZLHI60)s«TENP 

TEMP*SQRt <165.0t6/PPll/0.'8EO ~ 

ZLH(61)«21.CE0*rENP 

ZLm62>«-14.0E0*TEH» 

ZLN(63MrEKP 

TEHPaSQfT (11.0E)/FPI)/a. OEO 

ZLHf64l»63.0E'»*rEHP 

ZLMU5I«-70.0E0*TEHP 

ZtKI66)al5.0c.0*reNP 

ZLH(67)«ZLM(61I 

ZLH(&aiaZLH(F2l 

ZLH(&9)sZLHI63> 

TEPP*S(KT 11155. J’C/FPI f/2.0£0 
ZLM(70»«3.CE0*T£MP 
ZLMirifa-TEIIP 
ZLtH72l*-ZLMI54) 

?LM173>*-ZuH«55I 

2LM«74»«-ZL*U53> 

ZLm75l*-ZLM«56l 

ZLMI75 »rSOkT<J465.0rn/FPH/2.0E‘' 
7Lm77l»-ZLH!76* 


or!g;: ■ 

OF POOrl 


ZLMl7ft)«ZLH(49) 

ZL>M79|sZt.H(<»ei 

ZLHt60)sZLM(47) 

TEHP«S(MT 13003. SEO/l5t2.0EO*FPI> > 
ZLNfSl|s6.0E0*TEHP 
ZLH(82)«>2C«0EQ*TEHP 
ZLH(93|32lN(81) 

ZLHIo4)«SQRT(9009.0Ea/ 1128. 0E0*F FI ) > 

ZLH(85»«-ia.CEO»ZLM(94i 

ZLH(b6)s9.0E0*Zl.H(84) 

TENF«$QRT (619.0EC/(Z56.0EO*FPI)I 
ZLH(e7)>li.0E0»TEMP 
ZLHI88»«-66.CE0*TEMP 
ZLH(89t*ZlH(87) 

ZLHC90»«“TEMP 

ZLm91»«6.0E0*TEHP 

ZLNI92M-TEHP 

TEHP«SQRT a36S.qE0/U28.QE0»FPI) I 

ZLHI93|sll.0E0*rEHP 

7».M(94l«-33.bEb*TEMP 

ZL.M95»«9.0E0*T;MP 

ZLH(0&i«>3.0E0»rEMP 

TEHPs$Q9T(1369.0E07(512.0EO»FPI) ) 

ZLM(97)s33.0E0»TEMP 

ZcMI98»«-ZLM(97J 

ZLHt99J»-18.CEC*TEMP 

ZlMUOOIs*18.9EO*TENP 

ZLmiOi»*TERP 

ZLHI102»s-TEMP 

TEMP«SQRT<273.CE0/FPI )/8.0EC 

ZLH(lu3l»33.0E0*TEMP 

ZLMIl04|s-3b.bEb*TEHP 

ZLHa05»*5,0£0*TEMP 

TEHPsSQRT C13.0EO/FPI)/16.0EO 

ZLH(106|s231.0EO*TEHP 

ZLmi07l»-315.0E0*TEMP 

ZLH(10e)«109.0E0*TEHP 

ZLHCl09|i-5.0E0*TEHP ' 

ZLNIllOMZLNIiaS) 

ZLH(111I«ZLM(104) 

ZLM(112)«ZLNI109I 

TEXP*SOPI <13fc5,0E07(l28.OE0*FPI» » 

ZLMU13I=33.0E0*TEHP 

ZLMUl41*-18.3Eb»tEMP' 

ZlHI115»*TEmP 
ZLNI116» *-ZLM«94» 

ZLMU17l«-Zi.HC91» 

ZLMU16l*-ZtM(95» 

ZLM(ll9»«*ZtH(96l 
TE»1P»SO"T l819.0E0/FPi» /4.0EO 
ZLH(l20|«tl.0E0*TEHP 
ZLM(121Is>ZlHU20) 

ZLi«122>*-TEM? 

ZLMU23MTEHP 

ZLM(124|sZLH(86) 

ZL**ll25»*ZL‘«(85» 

Zu«<126»sZLM<84» 


ZLH (12 7MSQRM 3003.0 £0/(612. 0E0*FPII » 

ZUH(i 28 M- 15 . 0 E 4 *Zli(i 27 l 

ZLH(129)«-ZLH(128» 

ZLN(13(iM>ZlH(127) 

kETURN 

ENO 


SUSHOUTINE iXf AHOINS#LSfNSfLl#NT#NTPtNSP#LSPfNSPI„ 

C.« .GENERATE THE TRANSLATION INDUCED EXRANSION TERNS FOR THE EC> INTEGRALS OF 
C...THE FIRST RIND 
C 

IFLAG«0 

C 

LTNIN«NAXOfO»NT^LS-NS*lWl _ 

LTNAX«NINOaS-l#NT-ll41 
DO 50 LT«LTNIN»LTNAX 
NTN1N»NAX0C-LT«1#NS-2HS«LTI 
NTNAX>NINOfLT-l»NS*LTI 
HTNAX-NTNAX-NTNIN^i 

DO 40 NT*l»NTNAXfE 

C 

NXT«CLT^NT>NTNIN-2l/2 

NYT*(LT-NT-NTNIN)/2 

HZT-NT-LT 

C 

C _ _ _ 

LtFNIN-NAX6(0»NtFHSMSF-ll«l 
LTFNAX>NINOCLSF-l»NTF-im 
DO 30 LTPoLTFNINfLTFNAX 
HTFNIN-NAX0C-LTF^i#HSF-24LSFtLTFI 
NTPNAX-NlNOCLTF-XvNSF-LTFI 

NTPNAX-NTPNAX-NTPNIN^l 

DO 20 NTP«1»NTPNAX#2 
C 

NXTP«tLTP^HTP*NTPNIN-2l/2 

NYTP«CLTP-HTP-NTPNINI/2 

NZTP«NTP-LTP 

C . , 

c 

JFLAG-0 
NlNAX«24Ll-l 
DO 10 Ni«l#NlNAX 
C 

rENP»ECQCNXT#NYT»NZT#Ll-I#Nl-Ll#NXTP#NYTPfNZTPI 

IF (ABSCTENPI.lt. I.OE-IOI GO TO lo 
C 

IFLA6«IFLA6^1 
IF CIFLAG.GT.i) GO TO 5 
PRINT 520fU-l#NT-lfNTP-l 
5 CONTINUE _ „ 


4FLAG«JFLAG41 
IF (JFLAG.GT.il GO TO 7 

PRINT 530»NT-l»LT*lfNT4NTNlN^lfNXTfNYTfNZTfNTP-l#LTP-l» 
1 NTP4>NTPNIN-l»NXTPfNYTP9NZTP 

7 CONTINUE 
C 

PRINT SAOiNl-LlfTEHP 
C 

10 CONTINUE 
C 

20 CONTINUE 
30 CONTINUE 


40 CONTINUE __ _ _ 

50 CONTINUE 
C 

RETURN 

C 

520 F0RNAT(lH09l0Xf3HLl-f I292Xf3HNl«fI2»2Xf 3HN2«f 121 
530 FORNATC 15Xf5XfiH(f3I293H|/ (#3129 IHIf 32X93X9 1H(»3I293HI/( 9 312# IHII 
540 F0RNAT(42X92HN«9l293X9E15«BI 
END 


c::— - 

OF FOOk / 


OF pook q V 


SUBROUTINE INl(lNO(H$#LS» NS^NSRp LSRf HSRI 
C*..LOOP OVER THE TERNS IN THE ECP 1NTE6RHS OF THE FIRST KIND 
C 

IIHAX«NS4NSF*1 
DO 30 tl-ifUNAX 
C 

00 to NT«1»NS 

C 

DO 10 NTF-lfNSF 
C 

c 

C*«.6ENERATE THE TRANSLATION INDUCED EXPANSION TERNS 

CALL EXPAND(NSfLS#NS»LUNTfNTPfNSP#iSP»NSP| 

C 

10 CONTINUE 
20 CONTINUE 
30 CCNTINUE 
C 

_ RETURN 

C 

END 


OjO 


SUBROUTINE EXPANO(NS»LSf HSfLtUt UtNTtNT^tNSPtLS^tNSP) 
u...CENEMfe The TRiSsUffON InoucIC IxPiNStON TS»HS"T60 hT I5f 
C*..OF THE SECOND KIND 
C 


C 


C 


lFLtG >0 


CTNlNsNAXC 10 tNT»LS«NS*U *i 

LtNAX«NiNOaS-itNt*i)M 

00 70 LT>LTNIN,^rN«X 

NTNINaNAXO (*17 *l * HS-2* LS ♦Lf » 

NTHAX«NIN0<I.T«1|NS-ITI 

NTNAXaNTNAX-NTNlNti 

00 SO NT»l,HTNAi(,Z 


ohigi.'sal Pun-7 
PO.JS o’jiLrn 


NXI«(LT*HT*HTHIM*2I/2 

NTT»(tT-HT-HTNINI/2 

N2T«NT-LT 

C 

C 

LTPHrNiilAX0(6*NTP*L§P-f6P-lMl 
lTPHAXsNINO <lSP> 1»NTP-1I«1 
00 50 LTPaLTPHiMttTPNAX 
NTPNIN«NAXO »-LTP*l,NSP-2*LSPHTP) 
NrPNAX«HlNO(LrP«l»MSP»LTP) 
NTPMAXaMTPNAX-HTPHINtl 
00 AS HtP«itNfPNAXt2 
C 

NXTP* ILTP»HTP*Hf PHIN-21/2 

NTTP«iLTP-HTP-MTPHIN>/2 

NZTP»NTP-LTP 


JFwA6«0 

C 

H1MAX»2*L1-1 
00 3C Hi«l»HlHAX 
H2HAXs2*LC-l 

jO 2i M2*i.H2HAi( ' - -- - 

HMAX«2»L-1 
JO Vr H*1»1HAX 

r 

TE.il l«uCn(Ll-ltil-wl fNXT tNTT.NZT tL*ltH-U 
rEHP2*C;C0a2*l t«2-L2tNXTP,NrTPtN2TPtL-ltM-.l 

IF ICA 6 S ITEMPl) .LT.l, 5 c-id» .OR.IioSrtEMPai ,Lfii,irE-iDT~l 5 I 5 T 6 ' IB 

IFLA3«1FLAG*1 
IF tIFLAG.OT.n GO TO 5 
PRINT fcOJtt-lf Ll-i.w2-l»NT-l,NTP-l 
5 CONTINUE 


JFlAG«JFlAG+1 
aF UFLAG.GT.II go to 7 

PRINT 50?,NT-1 ,.T-l,NT»MTiIN-l iNXT f NTT ,N2T ti TP-1 1 LT» -1 , 

1 MTF*irPHii-i,NXT?tNtTP,Nzrp 

7 CONTINUE 

PPiNT 510»M1-Llf H-L,TE.iPl,T!MP2 

C . 

10 CONfINUE 
ec CONTINUE 
30 CONTINUE 

c 

c 

40 CONTINUE ... __ 

50 CONTINUE 
C 

c 

so CONTINUE 
TO CONTINUE 

C ..... 

RETURN 

C 

40 0 FOrP^T (IMO* 10Xo2HwS, 1 2 t2X 1 3NL1 «t 1 2 t 2X * 3ML2« # 1 2t 2* t3HNl»t I 2f 2Xt 
I 3MN2»oI2) 

5:C FO^MftTCl5X*5X#lM(t3l2f3Ml/Tf3l2tlHlf3?Xf3XolMl#312»?M>/lf3I2olMM 

51^ FO*^M4TnOX* 3MM1«, 12# 2X * 3MN2» t ^2» 2X ^2HH« 1 12f 3X^t£l5_» 0 S | 

ENO ■ 


OF POOR QUALITY 


SUBROUTINE INKlN0<NSfLSff1StNS»»f-SPfHSPj_ 

U.;.LbOP bitk tHE TERNS THE IfiP 1 We5r4LS OP TME SECONS KIsT 

c 

DO $0 
C 

L1HAK«U^NS-1 
00 40 Li«l*tlHAX 

c 

L2HAX«L^NSP-1 
00 3a L2«lfL2NAX 
C 

00 20 Nr«ltNS 

c 

30 ia NTP*i,MSP 

c 

CALL EXPAN0CNStLStN5fL«LitL2fNTtNrP«NSP»lSPtHSP> 

C 

10 CONTINUE 
20 CONTINUE 
30 CONTINUE 
40 CONTINUE 
50 CONTINUE 

RETURN 

C 

END 


